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ABSTRACT 

Computation of the terms of the vorticity equation 
is made at 300 mbs using the assumption that vertical mo- 
tions are negligible and approximating relative vorticity 
by geostrophic vorticity computed from finite differences. 
The validity of the computations is then tested, giving 
good qualitative agreement but poor statistical agree- 
ment. The isobaric divergence at 300 mbs as computed 
using 2 = Ve: V is compared with the value of 
divergence as obtained from the vorticity ren. Again 
there was only qualitative agreement. The magnitude of the 
lateral and longitudinal components of isobaric divergence 
were studied and it was found that the lateral component 
was dominant in 75 per cent of the cases studied. No rela- 
tionship was found between the longitudinal component and the 
total divergence. The synoptic aspects of the divergence 
fields at 300 mbs are discussed with respect to the #00-:85 
surface. 

The writer wishes to express his appreciation for the 


guidance and encouragement given him by Professor F. L. Martin 


of the U.S. Naval Postgraduate School in this investigation. 
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LIST OF SYMBOLS USED 
Rectangular (Cartesian) coordinates; standard 
coordinates with x axis toward the east, y axis 
toward the north, z axis toward the zenith. 


Rectangular (Cartesian) system of unit vectors 
along x, y and Z axes. 


Vertical component of absolute vorticity. 


= A three dimensional time rate of change operator. 
- Gradient operator on an isobaric surface. 
- Wind velocity vector. 


- Isobaric divergence 


E dp - Rate of change of pressure of an air 


particle with time. 


-Local rate of change operator with respect to 


pressure. 


- Horizontal component of wind vector. 


Iime rate of change operator along an isobaric 
surface. 


- Acceleration due to gravity (32.17 feet per second 


per second). | 


Coriolis parameter (2 Q.sin y ). 


Angular speed of the earth (7.292 x 107? radians 
per second). 


Latitude 
Height of an isobaric surface. 
Laplacian of height. 


Components of wind velocity along the x,y and z 
axes respectively. 


opace-averaged'height of an isobaric surface. 


Grid distance 2.5 of latitude -9.12 x 10? feet) 





df - Rate of change of coriolis parameter along 
E. y axis. 

Angle between wind velocity and the gradient 

of (mez). 

Density of the air. 

Natural system of unit vectors. 

Wind direction. 

Correlation coefficient. 

Covariance of stochastic variables x and y. 
standard deviation of a stochastic variable. 


Variance of a stochastic variable. 


Size of a statistical sample. 
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1. Introduction. 


Several studies of divergence in the atmosphere have 
been carried out by different methods and by different inves- 
tigators. Some of the most notable of these are Fleagle 12] ; 
Arnason [1] , and Landers [5 | . Their investigations have 
involved a three dimensional study of the distribution of 
divergence in the atmosphere. Landers [5] computed horizon- 
tal divergence using constant level charts and a "del com- 
puter" which consists of computing divergence from the wind 
vectors at the vertices of an equilateral triangle with an 
altitude of 249 of latitude. Arnason [1] computed the inte- 
grated divergence over a 12 hour period in isobaric layers 
with the use of vertical velocities, the local rate of 
change of thickness and the temperature advection.  Fleagle 
computed horizontal velocity divergence and mass divergence 
up to about 150 mbs and found a level of non-divergence at 
about 600 mbs. 

Ihis study is carried out at 300 mbs only and no attempt 
is made to extend it to three dimensions. The main objective 
of this study was to gain some practical knowledge of the 
behavior of divergence patterns and to apply the results to 
the prognoses of the 500-mb chart. Computations of the di- 
vergence and the individual rate of change of the absolute 
vorticity with time are made, and their relationship tested, 
Also a study of the relationship between the longitudinal 
component of the divergence and the isobaric divergence is 


made. 














2. Source of Data. 


À series of 300-mb fms, carefully analyzed by the 
staff of the U.S. Naval Postgraduate School, was used in 
making the computations. The analyses employed were those 
of 15002, 15 November, 0300Z and 15002, 16 November 1955. 
On the first chart there is a trough and a closed low 
at about 112° W and a ridge at about 849 W with consider- 
able meridional flow. At 0300Z on the 16th the trough and 
ridge have moved about 59 eagtward, retaining about the 
same intensity... By 1500Z on the 16th. the amplitude of the 
| trough-ridge system has decreased considerably and moved 
to 989 W and 729 W longitude so that in between the trough 
and ridge the flow is becoming zonal. Computations, which 
will be described in the following sections, were carried 
out between this trough and ridge. In addition some com- 
putations were made behind the trough on the 0300Z, 16 
November chart. A total of 138 computations of all five 


terms of Equation (7) below were made. 


3. Mathematical Analysis. 


The vorticity equation in isobaric coordinates is: 


× ت مل ۰ با .1.1 mu‏ )1( 


This is not in the form in which computations from a 





synoptic chart can be readily made. Since vertical veloc- 


ities are normally small, the second term on the right side 





of the equation may be ignored. Equation (1) then becomes 
the approximation to the vorticity equation as it is common- 


ly made in meteorology. 





i! DA ..5 (gv) 


Making the geostrophic approximation, 3:2 V! , and 
noting that the coriolis parameter does not change with 


time, Equation (2) becomes: 
(3) pae V2 +5)= - (Ev) v Y) 


In order to compute the individual rate of change of 
vorticity it would be necessary to construct isobaric tro- 
jectories of the air parcels. However, such computations 
may be avoided by making use of the fact thats 


DBF | 
5) De" TE AG. 


Eguation (3) then becomes: 
(5) (iv) + Ye (Vz) + vas : -( Se v) 


Since g is quasi-constant compared to the horizontal var- 
iation of V z , it may be brought outside of the differential 
operator in the second term of the left side of Equation (5). 
In order to get Vz ina form in which computations can be 


made, finite differences may be used to give: 





O Vs 








Combining Equations (5) and (6) and multiplying by fd^ gives: 


= Lou 
o ke etie (Pee = gen jen 


3 b c d € 


The equation is now in a form in which computations were 
made. The left hand side of this equation now represents 
6 አፈ and the right hand side represents -fd^$,(V, V). 
By neasuring wind velocities in feet per "T and grid 


distances in feet, the units of this equation can be shown 


to be feet per second. 


+, Computational Procedure. 


The quantity z may be readily computed using a grid 


as shown in Figure 1. 


Figure l 


Grid for computing space mean, 


e 


The value of 2d was chosen to be 59 of latitude. This 
value iz about half of the length of the short-wave trough 
so that in constructing a Z-chart, the short-wave troughs 


tend to be removed, but the long-vave troughs are unaf- 


fected. 





The construction of the z-chart where Z is zjt Zo*t 23+ 2 


—— 
(Figure 1), was accomplished by displacing two z-charts 10 
degrees of latitude laterally and adding graphically giving 
Z]* 23. Ihe two maps were then displaced 10 degrees of 
latitude northward and added giving z54z,. The 27+23 
contours are then added to the zo4zy contours giving hz. 
By dividing by 4 we then have a space mean map from which 
the actual map may be graphically subtracted giving iso- 
lines of (Z-z). The local rate of change of (z-z) was 
approximated by the average rate of change in the guantity 
over a l2-hour interval, as measured at any given station. 
The advection term | term (b)) was computed by measur- 
ing (1) the gradient of (Z-z) across selected stations 
over a distance of 59 of latitude, (2) the wind velocity 
in feet per second, and (3) the angle Y between the wind 
vector and the gradient vector. The magnitude of term (b) 
is then (VI Iv, (2-2)! \cos Yl. A diagram of the vectors 


involved is shown in Figure 2. 
v, (2 -2( 


۲ Y 


Figure 2 
Diagram of Angle ۷ 


All factors of term (c) with the exception of v are 


constant for a given latitude and can be computed for each 


2 





latitude circle. Taking the appropriate value of fa? B 
and multiplying by the northward component of the aii 
velocity gives the final term in the vorticity computation. 
Now by equation (7) the sum of terms (a), (b) and (c) rep- 
resent the quantity fd? DJ4 subject to the accuracy of the 
approximations made in EE cu vertical motions and in 
computing 3 (Z-z). 

Ne ከጨ have been devised for computing iso- 
baric divergencd | 3， 6| . It was decided to use the fol- 


lowing equation in the computation. 


. 3a ስ 
EXE V, Vs + 


cc 


A square grid 2% degrees of latitude on a side was used. 
Àn average u was computed on each side of the square and 
an average v at the top and bottom. As shown in Figure 3, 


the isobaric divergence is then E E 2T 
d a | 





Na) a] 


Figure 3 


Grid for computing divergence. 


Since the grid is sguare no spherical correction is 
necessary, and grids can be oriented along the major flow. 


If the values of isobaric divergence are then multiplied 
6 





by term (d) of Equation (7) the value, should be the same 

as the sum of terms (a), (b) and (c). The divergence was 
also computed by dividing the sum of terms (a), (b) and (oc) 
by term (d) of Equation (7). For the on in ۹5658 
paper this value of isobaric divergence will be referred 

to as effective divergence. 


Isolines of the computed values of fd“ DJ, are 
g Dt 


plotted in Figures hà, 8 and 12. The isolines of 
-fd^ JA (Y, Y) are drawn in Figures 5, 9 and 13. The field 
ی‎ divergence is plotted in Figures 6, 10 and 14, 
and isobaric divergence is plotted in Figures 7, 11 and 15. 
The labeled values of isolines of fa? D Ja, and 


g Dt 
-fay (Jp V) are times 1073 and the labeled values of the 
z pe 


fields of divergence are times 1077. In plotting the values 
of isobaric divergence in Figure 7 the values in four adja- 
cent squares were averaged and then plotted in the common 
point of the four squares.. This was done in an effort to 
smooth the divergence field but was not done in Figures 11 


and 15 because the averaging could mask the effects of short 


waves. 
5.. The Results of the Computations. 


The individual terms on the right hand side of Equation 
7 have order of magnitude - 107? feet per second. A typi- 
cal value of the change of (z-z) over a twelve hour period 


is 100 feet. When divided by 12 hours or 4.32 x 10 ‘seconds, 


N 


7 





<< (282) has a value of about 7 2.3 x Id feet per second. 


Ihe extreme computed value of 2- (Z-z) was 11.6 x 107 feet 
per second. Measuring A) over a distance of OG lati- 
tude (239 each side of the station) gives values of the order 
8 5.0 w 10”, With a wind of 120 feet per second and an 
angle Y whose cosine is 0.3, a typical value of "- (b) 

ES 1.8 x 1072 feet per second. The maximum magnitude of 

this term was 20.5 x 1072 feet per second. The guantity 

fa^ ranges from 11.4 x 1079 at 30° and 60° N to 


: at 1.59 N where it is maximum. When multiplied 


13.18 £ JO 
by a typical value of 100 feet per second for v, one obtains 
a value for term (c) of 1.1 x 1079 LOST. x 1073 feet per 
second. The maximum value of term (c) was 1.87 x en feet 
per second. Although term (c) is of the same order of magni- 
tude as terms (a) and (b), it is somewhat smaller in general 
and terms (a) and (b) give the largest contribution to the 
sum of the three terms of the left side of Equation (7). 

The value of f^d? is constant for a given latitude 
circle and ranges Teom 1.37 X 10* feet at 204 N to 3.66 x LO 
feet at 559 N, while (Z-z) ranges from -300 feet to 300 feet. 
Term (d) was positive in all the computations except one,.in 
which it was a small negative value, Of course all compu- 
tations were carried out in the broad scale flow between 
trough and ridge. Term (d) could very well be a negative 
value in a closed high pressure system where values of 
(z-z) might be as great as -hOO feet. 


Values of isobaric divergence were in general smaller 


8 





than those of effective divergence. However they agree 
fairly well qualitatively. The magnitude of the effective 
divergence agrees very well with that computed by Landers [5]. 
The average absolute value of the effective divergence on 

the three charts was 3.34, 3.70 and 1.22 all times 0 per 
second. A plot of divergence values versus height was made 
by Landers [5] and shows an increase of O.6 x 107? per 
second in the mean absolute value of horizontal divergence 
from 20,000 to 25,000 feet. Extrapolation of his curve to 
30,000 feet would give a value of 2.7 x 107? per second, and 
the mean absolute value ofeffective divergence of all the 
computation made in this study is 2.78 x 107? per second. 

The values of isobaric divergence have a maximum magnitude 

of about 3.5 x 107? per second while the mean absolute 

2 


value is about 1.0 x 10 ^ per second which agrees better 
with the values of divergence as computed by Arnason [1] š 
In the following section a statistical approach is made in 


comparing the two fields of divergence and the fields of 


faf Df, and -f82 4, (V, V) 
5 Dt E 


statistical Studies. 


Of the patterns of fd? D and -fd^ EV 
9 Det E pU = ) 
Figures 12 and 13 seem to give the best agreement. To test 
this agreement a correlation coefficient of the two quantities 
was found using values at 35 points in each figure. This 


correlation coefficient r is given by the folbwing equation 
9 





from Hoel [4]. 


(9) r =- cov (xy) 
ox Sy 


Cov(xy) is the covariance of x and y, the stochastic variables, 
and Sy and S, are the respective standard deviations.  Cov(xy) 
and the standard deviations were computed from the following 


equations given in [h]. 


(10) cov(xy) = E 2 xy- n Y j| 

(11) se [2-6 x) | 
2 

(12) ss = Liz (Ey) | 


The mean value of fa? DY was -2.06 x 103 feet per 
g Dt y 
second and the mean value of -fd2 fa (VV) was 0.1 x 10 3 


feet per second. The standard déviations were 6.38 x 107? 
feet per second and the 3,52 x በከት per second respect- 
ively. The correlation coefficient was found to be 0.31. 
With 90 per cent confidence limits the critical value of 
r is 0.29 so the correlation coefficient of 0.31 is signif- 
icant at this level of confidence. However, the value of 
r is not large enough to deduce more than that the two 
guantities are related. Some of the reasons for this lack 
of agreement will be discussed later. 

The same procedure as described in the previous para- 
graphs was used to compute a correlation coefficient between 


Figures lO and 11 which are patterns of effective divergence 


1O 








and isobaric divergence. It was believed that perhaps better 
agreement could be achieved between these two patterns since 
in the process of dividing Df, by -f some of the errors in 
using the geostrophic M mm might be minimized.  How- 
ever, the correlation coefficient between Figures 10 and 11 
was only 0.144 which is significant at 68 per cent confidence 
limits with a sample size of 52. 

The approximation of I (z-z) by a 12-hour average 
would tend to cause a computed effective divergence center 
to be east of its instantaneous position if the future 12- 
hour (Z-z) map is used and west of its instantaneous posi- 
tion if the previous 12-hour (Z-z) map is used. The plotted 
values of effective divergence and isobaric divergence were 
examined in an attempt to find some consistent displacement 
which would give a higher correlation coefficient. It was 
believed that the correlation coefficient could be improved 
by shifting the divergence center over Indiana in Figure 10 
into coincidence vith the divergence center over Michigan 
in Figure 11 while the convergence centers over lowa remained 
in coincidence. Actually, however, there was no improvement 
in the correlation coefficient. A value of 0.07 was com- 
puted for r with a sample size of 31 which is not significant 
at any reasonable level of confidence, 

Sutcliffe KA has stated that the vortex-tube term in. 
the vorticity equation tends to be offset by - T E v) 
so that a better approximation to the vorticity equation 


would be: 


E 





113) -É z-f(vev) 


À correlation coefficient was computed using this 
approximation and the computations from the 15002, 16 
Novenber chart, Using Equations (9), (10), (11) and (12) 
the correlation coefficient was computed to be -0.20h which 
is significant at 79 per cent confidence limits. The fact 
that the correlation coefficient is negative seems to indi- 
cate that ignoring f(9,:V )must overcompensate for the vortex- 
tube term. 

Using the eguation for isobaric divergence in natural 


coordinates: 


وم څه ۷ « ۷,۰۷ )14( 


it is seen that there are two components of the total di- 
vergence, a longitudinal component along a streamline and 
a lateral component normal to the streamline, It was 
also of interest to determine which of these components 
gives the greatest contribution to the total isobaric di- 
vergence. The rate of change of wind velocity along a 
streamline is most easily measured and was found at all the 
stations at which vorticity computations were made. This 
was done by measuring the change in wind speed along a 
streamline for a distance 5° of latitude across each sta- 
tion. If these values are compared with the effective 
divergence, the difference will be the effective lateral 
component. In most cases the lateral and longitudinal 


Mn 


— t 


1 








terms are of opposite sign. In 103 out of 138 computa- 
tions (approximately 75%) the lateral term was the domin- 
ant term. Therefore, most of the time the sign of the 
effective divergence is determined by the effective 
lateral component. 

It would be helpful if some relationship could be 
derived between effective divergence and the longitudinal 
component since the longitudinal component can be easily 


measured. The total divergence could then be determined 


۷ 


simply by measuring Em Using Equations (9), (10), (11) 
and (12) to determine a correlation coefficient between 
av and the effective divergence gave a coefficient of 


0.097 wnich is not significant at any reasonable level 
of confidence. Therefore, very little about the magnitude 
of the effective divergence can be deduced from the magni- 


tude of the longitudinal divergence. 


7. Sources of Error. 


Ihe discrepancy which exists between computed values of 


£a > and ESSEN) and between effective and iso- 
baric divergence is due partially to the fact that the vortex- 
tube term and the vertical advection of vorticity have been 
ignored. However, the discrepancy is greater than the values 
to be expected from the terms that have been ignored, Due to 
the lack of actual wind reports at 300 mbs, it was necessary 
to interpolate between the actual wind observations. Error 


could have been introduced here. Also a 109 error in the 





observed wind direction and a 5-knot error in wind speed 
may cause a maximum error of 3.0 x 10% per second in the 
computed isobaric divergence while the same error in vind 
direction and speed will cause only an error of 0.34 x noe? 
in the effective divergence. However, there is another 
approximation in the use of the vorticity equation. This 
was the determination of Se (2-2) as a l2-hour average 
using two consecutive (z-z)-charts. The magnitude of the 


errors due to this approximation is unknown. 


9. Synoptic Aspects. 


lt was difficult to determine the track of the diver- 
gence centers as many of them disappeared or were out of the 
area of computations in the 12-hour intervals. For example, 
in the l2-hour interval from Figure 6 to Figure 10, three 
divergence centers have appeared where there was previously 
only one, 

Ihe effective divergence patterns give good qualitative 
agreement with the synoptic situation. Fleagle [2] has found 
the 500-mb level to be approximately a level of non-divergence 
and the 300-mb surface is near the level of maximum divergence. 
Landers [5| also found the 300-mb level to be a level of maxi- 
mum divergence and the 500-nb level to be a level of least 
divergence. Arnason's results [1] did not show such a distri- 
bution of divergence, but it is believed that the isobaric 
divergence as computed at 300 mbs may give a good indication 
of rises and falls of the 500-mb surface. In Figure 6 the 

14 


areas of effective divergence between trough and ridge are 
relatively large and clear cut. With this divergence 
pattern aloft we would expect the trough to continue to 
move eastward, which actually occurs. Figure 7 does not 
give this same agreement as there is a convergence center 
centered over station 534 which is behind the ridge. This 
would indicate retrogression of the ridge which does not 
occur. lt may also be noted in Figure 6 that as the strong 
winds east of the ridge at 8h? W enter the cyclonically 
curved isohypses of weaker gradient there is a change from 
convergence to divergence. On the left hand side of Figure 
6 there appears to be a perturbation leaving the long-wave 
trough as indicated by the divergence center over Colorado. 
As yet there is no actual short-wave trough indicated at 
500 mbs. On the next 500-mb chart, 030042 16 November, a 
short-wave trough has appeared east of the long-wave trough. 
Figures 10 and 11 indicate closed divergence centers ahead 
of the short-wave trough and closed convergence centers 
behind it. With the short-wave trough in the broad scale 
flow,. the pattern of divergence is much more complicated and 
the gradients are much greater. On the West Coast in Figures 
10 and 11 there is a center of convergence located ahead of 
a short-wave trough. Twelve hours later the trough is prac- 
tically non-existant. The divergence center located over 
Oregon may be associated with the trough that appears off 
the West Coast on the 1500Z 16 November chart. In Figures 
14 and 15 the short-wave trough of Figures 10 and 11 has 

15 





moved into the ridge and weakened. Consequently the divergence 
and convergence centers that were associated with it are 
not as clear cut. The patterns have become elongated 
along the COS with considerable convergence east of the 
trough. Therefore, the trough may be expected to decrease in 
amplitude and the flow to become more zonal, which is verified 
on the 03007, 17 November 500-mb chart. 

Due to the fact that the surface maps are two and a 
half hours earlier than the upper-air charts, it is diffi- 
cult to compare divergence patterns with the surface maps. 
However,. in the few cases studied the low-pressure centers N 
at the surface move toward the centers of divergence at 
300 mbs. It is also interesting to note that areas of di- 


“vergence tend to overlay the fronts at the surface. 


9. Summary and Conclusions. 


This investigation of the isobaric divérgence patterns 
at 300 mbs has shown that the technique of computing diver- 
gence from the rate of change of absolute vorticity gives ' 
very good agreement in magnitude with the computations of 
isobaric divergence made by Landers [5] . The computations 
made in this investigation indicate that a mean absolute 
value of divergence at 300 mbs is about 2,8 x 7 per second. 
This value of effective divergence acting over a column of 
air 12,000 feet thick for a period of 12 hours will cause a 
fall in the 500-mb surface of about 200 feet. The values of 
effective divergence plotted in Figures 6, 10 and 14 give 


good indications of future areas of rise and fall at the 
16 





500-mb surface. 

Ihe computations of isobaric divergence using Equation 
(8) gave a smaller mean value, which gives good agreement 
with the magnitudes computed by Arnason [1] 。 Except in 
Figure 7, these values of divergence also give good indi- 
cations of future developments on the 500-mb surface. 


There were rather poor statistical correlation co- 


efficients between computed values-of fac D da and 
8 


2 
-fd” ŷ (y, v ) and also between effective and isobaric diver- 





gence. This poor agreement is believed due largely to errors 
in, and insufficient, wind reports at 300 mbs. Also the 


12-hour approximation of 2- (Z-z) introduces error, 

The lateral component of the effective divergence was 
the dominant component in 75% of the computations made. 
There was found to be no relationship between the effective 
divergence and either the lateral or longitudinal component 


of divergence. 
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